First Detection of CO in a Low Surface Brightness Galaxy 
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\0 . ABSTRACT 

We report on the first attempts at searching for CO in red low surface brightness 
galaxies, and the first detection of molecular gas in a low surface brightness (jUb(0) o & s > 
23 mag arcsec" 2 ) galaxy. Using the IRAM 30m telescope, CO(l-O) and CO(2-l) lines 
were searched for in four galaxies - P06-1, P05-5, C05-3, & C04-2. In three of the galaxies 
no CO was detected, to Tmb ~1.8mK (at the 3cr level). In the fourth galaxy, P06-1, 
both lines were detected. Comparing our findings with previous studies shows P06-1 
to have a molecular-to-atomic mass ratio considerably lower than is predicted using 
theoretical models based on high surface brightness galaxy studies. This indicates the 
N(H 2 )/ / T(CO)dv conversion factor for low surface brightness galaxies may currently 
be consistently underestimated by a factor of 3 - 20. 
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Subject headings: galaxies: individual: (P06-1, P05-5, C05-3, C04-2) - galaxies: ISM 
galaxies: evolution 



1 Based on observations carried out with the IRAM 30 m Telescope. IRAM is supported by INSU/CNSR (France), MPG (Germany) 
and IGN (Spain). 
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1. Introduction 

Recently it has become clear that there exists 
a large population of low surface brightness (LSB) 
galaxies which previously had escaped detection (e.g. 



[mpey ct.al. 1988; Schombcrt ct al. 1992; Morshidi 



Esslingcr, Davies, fc Smith 1999; O'Ncil, Bothun, & 



Cornell 1997a; |Q'Ncil, Hauscr, fc Bothun 2000| ). The 
importance of LSB galaxies in the local universe is 
emphasized by studies showing that the majority of 
the galaxies in the local universe are LSB galaxies. 
Realizing then that most galaxies are optically diffuse, 
it becomes important to understand these systems if 
we wish to understand galaxy formation and evolu- 
tion as a whole. This study is made difficult, though, 
by the lack of data on LSB galaxies. In particular, 
in spite of the importance of LSB galaxies there is at 
the present time little information about the molecu- 
lar ISM in these systems. 

To date four studies have been undertaken to look 
for CO in low surface brightness system s. The first 



two studies - Schombert et al. (1990) and dc Blok 



[fc van dcr Hulst (1998) - focused on galaxies with 



(B — V) < 0.6 and achieved 3a detection limits of 
T MB =8A - 24.3 mK (with the NOAO 12m telescope 
and the JCMT). The next two studies (Knezek 1993; 
Braine, Herpin, & Radford 2000) were of giant LSB 
systems. In all four cases no CO was detected. A 
fifth, very recent, study by Matthews (2000) is also 
well worthy of mention. This study was designed to 
look for CO in late-type, edge-on galaxies. CO was 
detected for a sample of galaxies, a subset of which, if 
de-projected from their edge-on orientation, may be 
LSB systems. 

The recent discovery of red (B — V > 1) LSB 
galaxies ( p'Ncil et al. 19971 ) has brought the idea of 
searching for molecular gas in LSB systems back to 
the forefront. As red colors should be indicative of 
reasonably high metallicities and therefore should in- 
sure against CO depletion (for example, as is seen in 
low metallicity dwarf galaxies), it seemed worthwhile 
to make a fifth attempt at CO detection in these dif- 
fuse systems. 

2. Observations 

Our targets were chosen from the catalogs of O 'Neil, 
Bothun, & Schombert (2000) (OBS hereafter). In 
an attempt to insure against any anti-correlation (or 
correlation) between the galaxies' H I and CO con- 
tent (i.e. Kenney & Young 1989) the total H I mass 



and mass-to-luminosity ratios of the chosen galaxies 
vary widely, from 1 - 75 x 10 8 M Q of HI and 1.3 - 
4 M0(HI)/L Q (B). Additionally, as mentioned above, 
to insure against CO depletion the four galaxies have 
(B-V) >1. 

We observed the CO(1-0) and (2-1) rotational 
transitions with the IRAM 30m telescope in the pe- 
riod from June 24 to 28, 2000. Table 1 lists the 
adopted positions (determined using the digitized 
Palomar sky survey plates) and heliocentric velocities 
(adopted from the measured 21-cm values of OBS) 
for our target sources. The beams were centered on 
the nuclei of each galaxy. Pointing was checked ev- 
ery hour and was found to be better than 5" ,and 
focus was checked every 3-4 hours. For each source 
both transitions were observed simultaneously with 
two receivers in each band. To each receiver we con- 
nected a 256 channel filter bank with channel width 
1 MHz, (670 kms -1 bandwidth, 2.6 kms -1 resolution 
at 3mm). For the source [OBC97] P06-1, which was 
known to have a wide HI line (OBS) we used two 
512 channel filter banks with 1 MHz wide channels, 
as well as two sections of the auto-correlator of width 
510 MHz and resolution 1.25 MHz. We observed us- 
ing the wobbling secondary with a maximal beam 
throw of 240". The image side band rejection ra- 
tios were measured to be > 30 dB for the 3 mm SIS 
receivers and > 12 dB for the 1.3 mm SIS receivers. 
The data were calibrated using the standard chop- 
per wheel technique ( Kutner fc Ulich 198l| ) and are 



reported in main beam brightness temperature Tmb- 
Typical system temperatures during the first two days 
of our observations were 200 and 450 K in the 3 mm 
and 1.3 mm band respectively. During the last two 
days the weather conditions deteriorated and no use- 
ful 1.3 mm data were obtained under these conditions. 

3. Results 

3.1. The Detection 

One of the galaxies observed, [OBC97] P06-1, had 
a clear detection with both the J(l-0) and J(2-l) lines, 
providing the first detection of molecular gas in a LSB 
galaxy (Figure 1). The observed main beam bright- 
ness temperature was 3.2±0.6 and 5.3±1.1 mK, af- 
ter smoothing to 40 kms -1 (for the J(l-0) and J(2- 
1) lines, respectively). This corresponds to an inte- 
grated flux density of 0.95 ± 0.22 and 1.1 ± 0.33 K 
kms -1 . (Errors for the flux density are computed 
from the rms error of the main beam brightness tem- 
perature and an assumed 40 kms -1 error for the ve- 
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locity width.) 

The IRAM J(l-O) beam-width is only 22", consid- 
erably smaller than both the 3.6' Arecibo beam used 
in the H I observations of OBS and the optical diam- 
eter (d 27 ) of 61.5" reported by O'Neil, et al. (1997a) 
(Table 2). It is therefore worthwhile to consider 
whether or not the IRAM beam has captured the ma- 
jority of the observable CO. As this is the first detec- 
tion of CO in a LSB system, we are forced to turn to 
studies mapping the CO distribution in HSB spiral 
galaxies. In 1989, using 2.6-mm CO emission obser- 
vations of more than 200 galaxies, Young & Knezek 
(1989) (hereafter YK89) concluded the average CO 
extents of spiral galaxies are typically one-half the 
optical radius. If P06-1 follows this trend, the J(l- 
0) beam only covers approximately 50% of the total 
surface area of CO while the J(2-l) line is reduced 
to measuring only 12% of the total CO emission. In 
agreement with this idea is the comparison of H I to 
CO velocity widths done by [Tutui fc Sofue (1999)| , 



which put our measured velocity width at 66% of the 
predicted value. As the rotation curves of LSB galax- 
ies can be found to be more slowly rising than their 
HSB counterparts (i.e. de Blok, McGaugh, & van 
der Hulst 1996; van Zee et al. 1997; [B waters, Madore 



& Trewhella 2000| ), this may be an indication that 



the beam has not fully encompassed the galaxy. It 
is therefore quite possible that we have only mea- 
sured roughly 1/2 of the total CO(l-O) and 1/10 of 
the CO(2-l) emission of P06-1. 

3.2. The Non-Detections 

In three of the four observed galaxies, [OBC97] 
C04-2, [OBC97] C05-3, & [OBC97] P05-5, no CO was 
detected, with 3<r non-detection limits ranging from 
T MS =1.5 - 3.9 mK (Figure 2 & Table 2). These 
detection limits are factor of 6 - 7 smaller than the 
rece nt work done by |de Blok fc van der Hulst (1998 



and Bchombert et al. (1990) . It is interesting to note, 
though, that the integrated intensity limits of our ob- 
servations are only slight ly smaller than in the ie Blok 
fc van der Hulst (1998) sample due to the narrower 
velocity widths of their sample. It should also be 
noted that the optical size of these galaxies, 36", 24", 
and 30"diameters for C04-2,C05-3, and P05-5, respec- 
tively, indicates that the majority of each galaxy was 
enclosed by the IRAM beam (for the J(l-0) observa- 
tions). 



4. H 2 in LSB Galaxies 

The observed integrated CO flux of the four galax- 
ies can be converted into an H 2 mass by assuming a 
value of N(H 2 )// T(CO)dv = 3.6x lQ w cmr 2 / {Kkms' 1 ) 
(sec below for a discussion on the conversion factor). 
Assuming the J (1-0) value is a more accurate mea- 
surement for P06-1 (as the beam size at 115GHz sam- 
ples a larger part of the galaxy), this gives M# 2 = 
11x10 s M Q , or M H2 /M HI = 0.2. Comparing our re- 
sults with those of other spiral galaxies indicates P06- 
1 to have a molecular-to-atomic gas ratio near the low 
end of the observed values (i.e. YK89). This is sur- 
prising as the red color of P06-1 (B — V = 1) typically 
would indicate a much higher value of M# 2 /Mhi (see 
below) . 

For the non-detections, the 33ct (J(0-1)) Mh 2 /Mhi 
limits range from 0.09 to 0.6. Comparing this with our 
measured values for P06-1 shows the non-detection 
limits to lie from 1/2-3 times the detected Mh 2 /Mhi 
value of P06-1. Taking into account the relatively 
high inclination of P06-1 (i=60°) compared with most 
of the other LSB galaxies which have been observed 
(«tj/p=33°) makes it plausible that only slightly more 
sensitive CO studies could considerably increase the 
number of CO detections in LSB systems. This ar- 
gument is supported further by the recent discovery 
of CO in a number of edge-on, late-type galaxies by 



Matthews (2000) 



Using the detected CO lines to estimate the H 2 
content of P06-1 and the other galaxies studied in 
our survey is a highly uncertain process. No molec- 
ular line is ideally suited for the determination of 
N(H 2 ), as it depends on knowing the CO «-> H 2 con- 
version factor X=N(H 2 )/ J T(CO)dv, having a galaxy 
which is fairly homogeneous in temperature and den- 
sity (at least for the CO molecule), and understand- 
ing the optical thickness of the system. Additionally, 
in low metallicity systems (as LSB galaxies appear 
to be), CO emission is a poor tracer of molecular 
mass (i.e. YK89; |Braine fc Combes 1992) Mihos, 



Spaans, & McGaugh 1999 - hereafter MSM). To allow 
for comparison between our findings and the work of 
Schombert et al. (1990) and de Blok et al. (19 98) we 
adopt the value given by |Sanders et al. (1986)| ofX 



N(H 2 )// T{CO)dv = 3.6x10^' cmT' A / (Kkms" 1 ) de- 
rived from observing high luminosity IRAS galaxies, 
albeit with the warning that our H 2 calculations are 
for comparison only, and could easily be off by a fac- 
tor of 10 (or more!) (i.e. de Blok & van der Hulst 
1998; MSM). 
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It should be noted that we are comparing the mea- 
sured values of CO, without accounting for the con- 
siderably smaller beam size used for the CO measure- 
ments (11" and 22") than for the HI measurements 
(3.6'). If the CO extent is 1/2 the optical radius of 
these galaxies, as is found by YK89, then our compari- 
son is accurate for the non-detections but low for P06- 
1. As this is the first detection of CO in a LSB galaxy, 
though, there is no way of determining whether the 
findings of YK89 hold for these systems. 

5. Discussion 

Assuming a Salpcter-like IMF, near-IR observa- 
tions by Bell ct al. (1999)| of red LSB galaxies in- 
dicate they should have metallicities in the ~Z© 
range. When combined with a theoretical study of 
the molecular ISM in LSB galaxies (MSM), this indi- 
cates the CO intensity of the red LSB galaxies should 



be above 1.0 K km s 



(It should be noted that the [Bell ct al. (1999 



(for t he J=l transiti on). 

pa- 



per accurately predicts the values found by YK89 
in their observations of HSB spiral galaxies.) The 
MSM paper also provides a prediction of the observed 
Mh 2 /Mhi for Z Q galaxies ranging from 1-3, de- 
pending on the structure of the galaxies' ISM. For the 
three non-detections in our sample, both the observed 
fT(CO)dv and Mh 2 /M H i fall well short of the pre- 
dicted value. Our one detection, though, matches the 
predicted J T(CO)dv value well but has a Mh 2 /Mjjj 
value a factor of 6 - 20 too small. 

A number of factors may be contributing to the 
discrepancy between our observe values (and upper 
limits) and those predicted by both MSM and Bell 
bt al. (1999)| : 

1. If the CO extent is considerably larger than the 
beam size, our observed values could be un- 
derestimated. It is doubtful, though, that this 
idea plays much of a role. Even if the observed 
/ T{CO)dv value of P06-1 is half the actual 
value, this could only increase Mh 2 /Mhi by a 
factor of 2, and would most likely increase it by 
considerably less. Additionally, the three non- 
detections have optical (d27) sizes of less than 
36", insuring that most of each galaxy would fit 
within the beam. 

2. If the metallicity of the LSB galaxies is co nsid- 
erably lower than that predicted by Bell et al. 



1999) , the galaxies will have fewer dust grains 



than average molecular content. Additionally, 
low metallicity results in a low UV attenua- 
tion. This implies that the H2 (which is strongly 
self-shielding) would not be much affected by 
the lower metallicity, while the CO would more 
readily be destroyed (i.e. Schombert et al. 1990). 
Thus X, the H2 to CO ratio, would increase. 
For our observed galaxies, the theoretical predic- 
tions of MSM show that decreasing the metal- 
licities to Z=0.3Z Q will both increase X by a 
factor of ~2 and decrease M# 2 /M#j by 5. Si- 
multaneously, fT(CO)dv is reduced to 0.1 - 
0.3 K kms -1 (depending on the galaxies' cen- 
tral surface brightness). With these values the 
upper limits on / T{CO)dv for C04-2 and C05- 
3 could readily fall onto the values predicted 
by MSM for M# 2 /M#j and could come con- 
siderably closer to their predicted J T(CO)dv 
values. The upper limits of P05-5, though, are 
still short of the predicted M# 2 /M#j value by 
a factor of 3. 

Attempting to match P06-1 into MSM's predic- 
tions through lowering its metallicity also has 
problems. As the predicted and observed val- 
ues of J T(CO)dv match, lowering the galaxy's 
metallicity to 0.3Z Q can allow the predicted and 
observed Mh 2 /Mhi values to match, but the 
observed J T(CO)dv value is then too high by 
a factor of 10. 

3. The third factor which should be considered is 
the homogeneity of the ISM. MSM show that 
the galaxies with 'clumpier' ISM have higher 
values of X than their homogeneous counter- 
parts, as the denser regions are better shielded 
from the background interstellar radiation field. 

4. The last item to consider in these calculations 
is the efficiency of the galaxies in converting 
atomic hydrogen into molecular hydrogen. If 
the velocity dispersion of LSB galaxies is higher 
than in their high surface brightness counter- 
parts, the formation and growth of molecular 
clouds could be adversely affected. This would 
result in lower molecular content in LSB sys- 
tems. This idea is supported by the fact that 
the one LSB galaxy which has been detected 
in CO also has the highest HI gas mass (and 
presumably the highest total mass) of any LSB 
galaxy in this study. 



to drive molecule formation, resulting in lower 
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As can be seen from the above arguments, no one 
parameter can readily account for both the observed 
jT(CO)dv values/limits and for the low M H JM HI 
values which are observed both in this project and 
in previous searches for CO in LSB galaxies. If 
MSM's prediction for J T(CO)dv is correct, then we 
are handed a dilemma - either the four galaxies ob- 
served in this project are homogeneous, in which case 
one of the galaxies (P06-1) has a metallicity 10 times 
that of the other three galaxies, in spite of their simi- 
lar colors, or the galaxies have a clumpy structure and 
their mctallicitics are similar and all are a factor of 3 



- 1 ti mes fainter than that predicted by [Bcll| ct al. 
(19§9| In the latter case the models of MSM pre- 
dict that P06-1 has considerably more high density 
regions than the other three homogeneous galaxies. 

Assuming the MSM values for / T(CO)d v are cor- 
rect, at least to within a factor of ~3 (to account for 
possible variations in the HI^H 2 conversion rate), we 
have to face the issue that the previously used value 
of N(H 2 )/ / T(CO)dv is low by a factor of 3 - 20. due 
to the large number of uncertainties in X, it is quite 
likely that this interpretation is correct, at least to 
some degree. It should be emphasized, though, that 
such a high value of X may be indicative of a different 
evolutionary path or processes in LSB galaxies. 

Clearly there are still a large number of questions 
yet to be answered about these systems before a true 
understanding of their stellar structure and evolu- 
tion can be found. Perhaps the most important is- 
sue is a determination of the metallicities of the red 
LSB galaxies so that more accurate comparisons be- 
tween these objects and the models of MSM can be 
obtained. Regardless, though, two facts seem to be 
clear from our CO observations. First, the observed 
CO flux intensities of LSB galaxies are well below pre- 
dictions made which assume LSB and HSB galaxies 
to evolve in an identical fashion. Second, even the 
higher CO/H2 conversion factor used by Schombert 
et al. (1990) and de Blok et al. (1998) is most likely 
too low by a factor of at least two, resulting in higher 
than previously believed limits on M H2 in LSB galax- 
ies. 

PH acknowledges partial support from NSF-EPSCoR 
grant EPS-9874782. We also thank the anonymous 
referee for suggestions which improved the manuscript. 
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Table 1 
Observed Sources 



Source 


a(J2000) 


<5(J2OO0) 


VHel 


Mb(0) 




B-V 


f20 


M H i 


i 


d27 




h m s 


o / // 


kms -1 






km s~ 1 


x 1O 8 M 


o 


// 


C04-2 


08 23 29.3 


21 36 45 


5168 


24.0 


24.2 


1.3 


103 


4.2 


31 


35 


C05-3 


08 27 44.4 


20 28 50 


12940 


23.9 


24.0 


1.2 


130 


26. 


20 


24 


P05-5 


23 19 34.8 


08 45 41 


3177 


24.4 


24.6 


1.2 


71 


0.6 


37 


28 


P06-1 


23 23 32.6 


08 37 25 


10882 


23.2 


24.2 


0.94 


430 


74. 


66 


62 



Information taken from OBS, O'Ncil, Bothun, & Cornell (1997a), and O'Neil et al. (1997b) 



Table 2 
Line Parameters 



Source 


Transition 
CO 


(mK) 


/ T M b<1v 
(K km s- 1 ) 


M H2 
xlO 8 M 


M H 2 

M H i 


VHel 

(km s -1 ) 


Width 

(km s- 1 ) 


C04-2 


1-0 


< 0.5 a 


< 0.096 b 


< 2.3 C 


< 0.55 




(103) 


C05-3 


1-0 


< 0.6 a 


< 0.14 b 


< 2.2 C 


< 0.086 




(130) 


P05-5 


1-0 


< 0.7 a 


< 0.12 b 


< 0.11 c 


< 0.18 




(71) 




2-1 


< 1.3 a 


< 0.21 b 


< 0.047 c 


< 0.078 




(71) 


P06-1 


1-0 


3.2 ± 0.6 d 


0.95 


11. c 


0.15 


10904 e 


302 f 




2-1 


5.3 ± 1.2 d 


1.14 


3.2 C 


0.043 


10903 c 


216 f ' 



a For non-detections, Tmb is the la r.m.s. noise smoothed to a resolution of 40kms . 
b Non-dctcction limits are Ico = 3Tmb^= (N = the number of channels). 

<=M H2 = 5.82 [(tt/4) dflco] (Assumes N(H 2 )/ J T(CO)dv = 3.6 X 10 20 cm- 2 / (Kkms- 1 ).) 

d T_MS wa s computed as J Tmb^v divided by width. 
° Velocity of central line channel. 

f Width is given by number of channels above 3 a noise level. 
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0.005 - 



CO(1 -0) 



0.005 
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C0(2-1) 



0.005 - 



0.005 
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1.05 10 1.1 10 1.15 10 

Vi ( krT1 s_1 ) 



Fig. 1.— IRAM 30m CO(l-O) and (2-1) spectra 
for the low surface brightness galaxy P06-1. The 
data have been smoothed to a resolution of about 
40kms~ 1 . The horizontal bar indicates the extent of 
HI as measured by O'Neil et al. (2000). 




-200 



200 



Vi ( km s 1 ) 



Fig. 2.— IRAM 30m CO(l-O) spectra for three low 
surface brightness galaxies where we derived upper 
limits on the CO (1-0) line strength. The data have 
been smoothed to a resolution of about 40 km s" 1 and 
are plotted as offsets from their heliocentric velocities 
listed in Table 1. The horizontal bar indicates the 
extent of HI as measured by O'Neil et al. (2000). 
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